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We characterize the long range dipolar scattering in 2-dimensions. We use the analytic zero en- 
ergy wavefunction including the dipolar interaction; this solution yields universal dipolar scattering 
properties in the threshold regime. We also study the semi-classical dipolar scattering and find 
universal dipolar scattering for this energy regime. For both energy regimes, we discuss the validity 
of the universality and give physical examples of the scattering. 
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Many theoretical proposals are based on dipolar gases 
in 2-dimensional (2D) geometries. For example, there are 
predictions of Wigner crystals [l[ , anisotropic solitons , 
and phonon instabilities Such theories show dipolar 
systems will lead to exotic and highly correlated quan- 
tum systems. Reduced dimensionality also offers another 
level of control to exert over ultracold matter. In pursuit 
of such quantum systems, there has been exciting ex- 
perimental progress in the production of polar molecules 
[J, l5( and in the production of quasi-2 dimensional ultra- 
cold gases [f| 0] ■ This makes it seem that the experimen- 
tal production of 2D dipolar gases is at hand. However, 
there is no simple understanding of dipolar scattering in 
2D, even for the case of the long range scattering. Such 
an understanding is an important first step in the study of 
these quantum systems. In this paper we study the prop- 
erties of long range dipolar scattering in 2D and present 
estimates of the scattering cross section. We consider the 
scenario when the polarization of the molecules [d = z) 
is perpendicular to the plane of motion (p — \J x 2 + y 2 , 
d-p — 0). We do not consider the complications of trans- 
verse confinement. Under these simplifications the dipo- 
lar interaction is: Vdd = d 2 / p 3 , where d is the magnitude 
of the induced dipole moment. The anisotropy of the in- 
teraction has been removed and the interaction is purely 
repulsive. 

To understand the scattering we use the dipolar length 
scale: D = pd 2 /h 2 where p is the reduced mass. In po- 
lar molecules this length scale can be quite large, or- 
ders of magnitude larger than the range of the short 
range interaction, po- Using D to rescale the 2D radial 
Schrodinger equation and expanding in partial waves, 
V>(p, <p) = E m e lm *<b m (p)/^p, the result is: 

/ d 2 m 2 - 1/4 2 \ , , N .. , . 

( -^2 + pT~ + ^3 ) = (Dk) 2 cl> m (p), (1) 

where p = p/D, k 2 = 2pE/h 2 , and E is the scattering 
energy. If p = po/D <C 1 then the only degree of free- 
dom is Dk, so once Eq. |T]) is solved for a given Dk the 
resulting scattering will apply to any quantum mechani- 
cal 2D dipolar system; this is universal dipolar scattering 



We have analyzed Eq. ([I]) and characterized the long 
range 2D dipolar scattering. We present analytic esti- 
mates of the scattering in both the threshold and semi- 
classical limits. In contrast to 3D, there is a diagonal 
s-wave dipolar interaction, and this leads to universal 
dipolar scattering in the threshold regime. We refer to 
the isotropic m = partial wave as s-wave. For thresh- 
old scattering in 3D, it is required that a s <SC D to have 
universal dipolar scattering U3] ■ This is challenging 
because a s , the s-wave scattering length, depends sensi- 
tively on the details of the short range interaction. How- 
ever in 2D, the repulsive dipolar interaction prevents ul- 
tracold particles from reaching the short range interac- 
tion where it can acquire a non-universal phase shift. To 
have universal dipolar scattering we require only that the 
range of the short range interaction is much smaller than 
the dipolar length scale (po <C 1). This criteria does not 
rely on serendipity to be met. 

In this regime, the s-wave scattering properties are 
derived from an analytic zero energy solution to the 
Schrodinger equation including the dipolar interaction. 
This solution shows that dipolar interaction results in an 
s-wave scattering length of 6.3441?. We also present the 
Born approximation which predicts the threshold scat- 
tering for non-zero partial waves. In the large Dk limit, 
the semi- classical universal dipolar scattering cross sec- 
tion can be estimated by the Eikonal approximation and 
scales as y/ D jk. These estimates offer accurate predic- 
tions of the scattering in their respective energy regimes. 

To monitor the scattering we use the scattering cross 
section, which in 2D for distinguishable particles is a — 
!Er=-co sin2 (^)' where 5 m is the scattering phase 
shift for the m th partial wave [TJ [H, [l3| . For identi- 
cal bosons (fermions) the sum is restricted to even (odd) 
m and there is an additional factor of 2. Note that in 
2D, a has units of length. We numerically solve Eq. ([1]) 
and obtain the scattering phase shift by matching <p m 
to the free solution at large p: ^/kp[cos(5 m )Jrn(kp) — 
sm.(8 m )N m (kp)], where J rn /N m are regular Bessel func- 
tions of order m. An interesting fact of 2D scattering is 
that the s-wave, cross section diverges as k — > 0. To illus- 
trate this, we solve the s-wave scattering with no dipolar 
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FIG. 1: (Color Online) (a) The 2D scattering cross section 
is shown for bosonic KLi (red squares, at), fermionic RbK 
(blue circles, erf), and distinguishable RbK-KLi (black +, ad) 
for D — 5.29 x 10 _6 cm. (b) The dimensionless ka is shown 
with three analytic estimates; they are the Born Approxima- 
tion (blue solid line), semi-classical (solid green line), and the 
analytic s-wave threshold scattering. The symbols are data 
from (a) with their energy range extended to show the low 
energy trend of the scattering. 



interaction and a hard wall at po- This yields 

<Jhw{po) = y sm 2 (S hw ), 
k 

twi(bhw) = J Q (kp )/N (kp ). 



(2) 



As k goes to zero the phase shift becomes logarithmically 
small because No oc ln(fcpo); however this is not quick 
enough to counter the diverging 1/k. 

We have plotted the cross section for three physical sys- 
tems. In Fig. [Ha) we show the cross section as a function 
of energy for bosonic LiK (a^), fermionic RbK (cr/), and 
distinguishable scattering between LiK- RbK (ad)- We 
have used D = 5.29 x 10 _6 cm (10 3 ao); this value of D can 
easily be achieved for each system and is discussed be- 
low. Note the divergence of the cross sections containing 
s-wave contributions, and the Wigner suppression of cr/. 
In fig. [TJb) we show the dimensionless quantity ka as a 
function of Dk. ka is proportional to the scattering rate, 
and (Dk) 2 is twice the ratio of the kinetic energy and the 
dipolar energy of the system, 2E / (d 2 / D 3 ) . The symbols 
are the data from (a), but the energy range has been ex- 
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FIG. 2: (Color Online) The analytic s-wave zero energy wave- 
function (solid black line) and the zero energy free solution, 
(red squares) are shown. Also plotted is a wavefunction in the 
threshold regime (Dk = 5 x 10" 3 , blue x), which agrees well 
with <j>o a t small p/D. The inset shows these wavefunctions 
close to the origin. 



tended to show the low energy trend of the scattering. 
In fig. [lib) we also show three analytic methods to es- 
timate the scattering cross section: Born Approximation 
(solid blue line), threshold s-wave scattering (black trian- 
gles), and semi-classical (solid green line); each of these 
estimates is detailed below. 

To obtain the threshold scattering properties of the 2D 
dipolar problem, we use an analytic zero energy (k = 0) 
solution to Eq. (pj: 



(3) 



where Ki m /li m are modified Bessel functions. C is de- 
termined by the boundary condition: 4>q(po) = 0- Con- 
sidering the s-wave case in the limit of small po, C be- 
comes small and is proportional to e _2 v' 8 /^ . In this limit 
4>o(p) = 2VpK2m(y/8/p) and this wavefunction is plot- 
ted in Fig. [2] as a solid black line. We have also plotted 
<t>Q (blue x) for Dk = 5x 10~ 3 , and it agrees well with 
4>q at small p/D, as is shown by the inset. 

To exploit the analytic solution, we find the form of the 
free solution, <fi Q , which is analogues to the zero energy 3D 
wavefunction: ripQ D (r) — r — a s , which has the property: 
ipQ (a s ) — 0. In 2D the scattering length is defined as 
where zero energy free solution is zero: cf>l = \fphi(p/a s ) 
[l4| . To extract the s-wave scattering length, we take the 
large p limit of cAg and find: 



fp [-2 7 - ln(2) + ln(p)] = ^fphx(p/a 3 ) 



2 7 +ln(2) 



D 



6.344. 



(4) 



where 7 ~ 0.577... is Euler's constant. Note 4>q in the 
large p limit is the zero energy free solution, and the free 
solution is plotted in Fig. [2] as red squares. The inset 
shows the wavefunctions near the origin. The thresh- 
old cross section for 2D dipolar scattering can be ap- 
proximated by the cross section of a hard wall at a s , i.e. 
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(5) 



This is a remarkable result in both its simplicity and ac- 
curacy. In Fig. [3]^a) we have plotted k<7d (black +), fc<r& 
(red squares), kahw (black triangle) for both the bosons 
and distinguishable case, and kcrsc (solid green line) as 
a function of Dk. From this figure we see Eq. ([5]) repro- 
duces the scattering well in the threshold regime when 
Dk < 10~ 2 . 

To predict the threshold scattering of non-zero partial 
waves we use the Born Approximation (BA) [ll|, [ijj , and 
find 



{Dk) 2 



k (to 2 - 1/4) 2 



(6) 



This offers a very good estimate of the non-zero partial 
waves scattering cross sections. For identical fermions, 
the BA gives kaf A ~ 2%Dk) 2 . We have plotted the BA 
(blue solid line) in Fig. [IJb) for k<7f (blue circles). In 
the small Dk limit, the BA reproduces the kaj well. In 
Fig. [T](a), w e see that the fermionic cross section (aj) is 
Wigner suppressed, as expected This is in contrast 
to 3D where all non-zero partial wave become energy 
independent as energy goes to zero [Icj . The analytic 
solution presented by Ref. [l(| produces the results in 
Eq. © when applied to 2D. 

To estimate the transition out of the threshold regime 
for to = 1, we can also equate the centrifugal and dipo- 
lar terms and find the energy of that length scale to be 
roughly Dk ~ 0.3. This is a reasonable estimate of when 
the scattering leaves the threshold regime, and is seen 
in Fig. [T] (b). The transition from the threshold regime 
does require higher energy collisions. D can be larger, 
which activates more partial waves. 

The BA breaks down when Dk becomes large and 
we enter into the semi-classical regime. In this limit 
the scattering is no longer dominated by a single par- 
tial wave, and to understand the scattering in the 
semi-classical limit one must obtain the total cross sec- 
tion. To do this, we estimate the scattering ampli- 
tude with the 2D Eikonal approximation [l7j : = 
IT 00 dbcos(kb9)[e lS " k sm(5 e ik)], where the phase is 



-ik2 
Seik 



Ujp) 



, and b is the impact param- 

eter. For the dipolar interaction U(p) = 2D/p 3 , the 
phase shift is —Dk/(bk) 2 . Using the optical theorem, 
fj = \/8Tr/kIm[f(0)] we find the total semi-classical 
cross section to be 



(Jsc = tV nDk. 
k 



(7) 



To illustrate the universal behavior of the scattering, in 
Fig. [21(b) we plot ko\, (red squares), kud (black +), ka f 
(blue circles), and analytic estimate from Eq. ([7]): kasc 
(green solid line) as a function of Dk. At large Dk, ka 
clearly shows the semi-classical V Dk scaling, kab and 




FIG. 3: (Color Online) (a) The low energy s-wave scatter- 
ing for both kod (black +) and kaf, (red squares) with kasc 
(solid green line), and the kahw{a s ) (black triangle) for both 
the bosonic and distinguishable cases, (b) The high energy 
behavior of kof (blue circles), kab (black +), and kod (red 
squares) with kasc (solid green line). The V Dk scaling of ka 
is clearly shown. 



kaf oscillate out of phase about ka^. Fig. [3] (b) shows 
that the semi-classical approximation offers a very good 
estimate of the scattering when Dk > 1. 

We have characterized the long range scattering as a 
function of Dk in the po « 1 limit. This ensured that 
deviations from universal behavior due to the short range 
interaction were negligible. But this leads to an impor- 
tant question: when can the scattering be impacted by 
the short range interaction? We have surveyed this ques- 
tion and have found that for po less than 0.1 the short 
range plays little to no role in the scattering, even for 
Dk <C 1. When po exceeds this value, then the short 
range can resonantly impact the scattering and spoil the 
universal threshold scattering. Once po > 1, the scatter- 
ing cannot be characterized by dipolar scattering alone, 
and knowledge of the short range is required. A more 
thorough treatment of the short range impact on 2D 
dipolar scattering will be presented in the future. 

An aim of this work is to facilitate estimates of the 
cross section for 2D dipolar scattering. This can be done 
with Fig. |4j In (a) we have plotted k" 1 as a function of 
energy for for LiH (black solid line), KLi (red squares), 
RbK (green dashed line), and RbCs (blue circles). With 
this and the value of D one can readily determine the 
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FIG. 4: (Color Online) (a) k' 1 is shown as a function of en- 
ergy for several 1 E molecules: LiH (black solid line) , KLi (red 
circles), RbK (green dashed line) and RbCs (blue squares), 
(b) D/Do is shown as a function of electric field. With these 
two curves one is able to estimate the scattering cross section 
for physically reasonable energies and electric field values. 



scattering character and estimate a value of the scatter- 
ing cross section. 

To estimate D, we use Fig. H](b) where we have plotted 
D I Dq as a function of electric field. Dq is the length scale 
obtained by using the bare dipole moment, d 2 m/h 2 . For 
the 1 T, molecules [1, the polarization of the molecules 
depends on the strength of the external electric field (8) 
in relation to the rotation energy of the molecules; there- 
fore the field can be characterized by B/d^. Fig. 2] (b) 
shows D I Dq which depends on the induced dipole mo- 
ment squared (d 2 ) as a function of electric field over B/d^. 
Note that the D is only a fraction of this length. A few 
examples are: 

RbCs RbK KLi LiH LiCs 

Do/10- 5 (cm) 28 3 45 20 317 
£ (kV/cm) 0.8 3.9 4.2 74 0.5 
Many others are given in Ref. U For mixtures, 



such as LiK-RbK, we need to determine the induced 
dipole moment of each molecule at a particular field: 
D = did2m r /h 2 . To offer a concrete example, con- 
sider fermionic RbK at lfiK, where fc -1 ~5x 10~ 6 cm. 
In Fig. Q] (a) we used D — 5.29 x 10~ 6 , this requires 
E — 1.25£ = 4.9kV/cm. Then Dk ~ 1, and we use the 
semi-classical estimate and find ka ~ Ay/irDk ~ 7. This 
estimate yields a ~ 3.5 x 10~ 5 cm, which is quiet close to 
the numerical value in Fig. [Ha). 

In conclusion, we have characterized the long range 
2D dipolar scattering and have offered several means to 
estimate the scattering cross section. In the universal 
threshold regime we have found an analytic solution to 
the interacting 2D radial Schrodinger equation. From 
this solution, we have extracted the scattering length for 
2D dipolar scattering: a s = 6.344D. Additionally we 
have explored the semi-classical universal scattering, and 
th e res ults are that the scattering cross section scales as 
y/D/k, from Eq. Q. This regime can be reached by 
not only increasing temperature, but also by increasing 
D. The strong dipolar interaction activates many par- 
tial waves to scattering and the scattering can be semi- 
classical even at ultracold temperatures. An important 
application of this work will be the use of a s in many body 
theories to correctly account for the dipolar interaction 
[3 [3 H3| • Future work will explore the scenario where 
d ■ p 7^ and the inclusion of transverse confinement in 
the scattering. 
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